Columnar defects acting as passive internal field detectors 



O 

o 

(N 
o 

o 

o 

m 



A.V. Silhaneki*, L. Civale\ and M.A. Avila^t 
^ Comision Nacional de Energia Atomica- Centra Atomico Bariloche and Instituto Balseiro, 8400 Bariloche, Argentina. 
^Instituto de Fiaica "Gleb Wataghm", UNICAMP, 13083-970, Campmas - SP, Brazil. 

(February 1, 2008) 

We have studied the angular dependence of the irreversible magnetization of several YBa2Cu307 
and 2H-NbSe2 single crystals with columnar defects tilted off the c-axis. At high magnetic fields, 
the irreversible magnetization Mi{QH) exhibits a well known maximum when the applied field is 
parallel to the tracks. As the field is decreased below H ~ 0.02_ffc2, the peak shifts away from the 
tracks' direction toward either the c-axis or the ab-planes. We demonstrate that this shift results 
from the misalignment between the external and internal field directions due to the competition 
between anisotropy and geometry effects. 
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I. INTRODUCTION 

It is a well established fact that the presence of colum- 
nar defects (CD) in high temperature superconductors 
(HTSC) enhance the critical current (Jc) due to the 
strong pinning and the inhibition of thermal wandering 
when the flux lines lay into these tracks .Elia The direc- 
tional pinning produced by these correlated structures 
becomes.^ evident when the angular dependence of Jc is 
studied.l3'Q In the last decade several works have shown 
that a sharp peak in Jc(0) appears when the applied 
field H is |alisj:iejd with the direction (6c) of these lin- 
ear defectil'BETEl (Here 9 and Qd are the angles formed 
by the crystallographic c-axis with H and the CD respec- 
tively). However, this behavior only holds when H is high 
enough to ensure that the average vortex direction is par- 
allel to H. At lower fields, both material anisotropy and 
geometry effiects become relevant and modify the vortex 
orientationjj which consequently may not coincide with 
that of H. Since maximum pinning occurs when the vor- 
tex orientation (given by the internal field B rather than 
by H) is aligned with the CD, Jc(Q) should maximize at 
an angle Qmax ^ ©d- In other words, any misalignment 
between B and H manifests itself as a shift in the angular 
position of the peak in Jc(©) with respect to the tracks 
direction. 

In fact, we have recently shown in YBa2Cu307(Y:123) 
that, if the anisotropy effect dominates over the geome- 
try effect (as occurs in most cases for HTSC compounds, 
as we will show below) for fields lower than \QkOe the 
peak progressively departs from the tracks orientation 
and shifts toward the c-axisci (|6maa;| < |0d|)- On 
the contrary, if the geometry effect dominates over the 
anisotropy effect (as would be the case in more isotropic 
materials), the maximum in Jc(0) should move toward 
the ab-plane (|6maa;| > |0d|)- Although the latter ^ 
feet has not been seen until now, Candia and CivaleE2l 
have experimentally demonstrated that in the isotropic 
Pbo,gTlo,i alloy, the sample shape determines the inter- 
nal field direction. In that study it was shown that, re- 



gardless of the direction of H, at low fields the flux lines 
remain almost locked to the sample normal because the 
system gains energy by shortening the vortex length. 

In this paper we present a study of the angular de- 
pendence of the irreversible magnetization in several 
YBa2Cu307 and 2H-NbSe2 single crystals with a single 
set of CD tilted off the c-axis. In particular, we study 
the misalignment between the external and the inter- 
nal field directions for several sample aspect ratios and 
anisotropics. We use the uniaxial pinning of the CD as a 
detector of the vortex orientation in the bulk of the sam- 
ples. We clearly demonstrate the influence of the crys- 
tal shape in the determination of the vortex direction 
and find that the angular behavior of the critical current 
is well described by the competition between material 
anisotropy and sample geometry. 



II. THEORETICAL BACKGROUND 

A complete description of the vortex lattice be- 
havior in a superconducting material should include 
geometry effects, mass anisotropy, vortex-vortex and 
vortex-defects interactions. Whenever the system is 
in thermodynamic equilibrium, the internal field 
is determined by minimization of the free energ; 
G(B) = F(B) - ^ + (MzH)M. In this expression the 
magnetization M and the applied field H are related by 
H = B — 47r(l — j>)M, where i> is the tensor of demag- 
netization factor. The components of i> at the sample 
principal axes are {iyx,i^y,Vz), with + i^y + = 1- 
We adopt the notation that z coincides with the crys- 
tallographic c-axis, and that the x axis is perpendicular 
to both c and H. Standard minimization of G(B) with 
respect to By and gives, 

dF B, Bi-Hi 

■ ~-y,z. (1) 
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where 



In the intermediate fields regiwie {Hd -C -ff ^ Hc2) 
the free energy i^(B) is given bj^. 
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where 0b is the direction of the internal field, cq^ = 

\/ cos^ 8b + sin^ b , the anisotropy e — mab/mnc < 1, 
and \ab is the penetration depth for H ||c. 
By replacing (||) into (|^) we obtain 

■ /p, /K,t^z,e) sin(29B) lnfe + 1 

where flvy, i^z, e) = {1 — v^) — {i — Vy)^^ and the reduced 
field = H/Hc2{Qb,T). 

The result (||) only assumes uniaxial anisotropy and 
the coincidence of one principal axis with the c-axis, and 
it shows that under those very general conditions the mis- 
alignments due to both mass anisotropy and sample ge- 
ometry have the same field and temperature dependence. 
The function /(i>, e), which contains the combined effects 
of geometry and anisotropy, is the key ingredient of the 
low field behavior, as its sign determines whether Qb 
leads or lags behind 8. 

To be more specific, let's consider the typical plate- 
like shape of all the single crystals of both Y:123 and 
NbSe2 used in this study, with thickness t along the 
c-axis much smaller than the lateral dimensions 
and L 

Vy = t/Ly, thus i>x, Vy, (1 — Vz) ^ 1 
strongly anisotropic and the crystal is not too thin, then 



y. To a first approximation = t/L^ and 

If the material is 



(1 



> (1 



thus / > and Ob > 9. We 



will call this the "anisotropy-dominated" situation. In 
contrast, for thin enough samples of a not too anisotropic 
material {1 — Vz) < (1 — ^'y)e^, so Qb < ©• This is 
what we will call the "geometry-dominated" case. The 
extreme limit of this case, with an infinite slab {i>x — 
Vy = 0) and ignoring the anisotropy, has been discussed 
by Klein et al.cl. It is also worth to note that for an in- 
finite cilinder with axis perpendicular to H, where the 
geometry effects are expected to cancel out, i/^, = and 
t'l/ = J^z = 5, thus / oc (1 — e^) and^.(^) reduces to the 
well known expression for the bulk.Bu 

Eq. allows us to determine which should be the vor- 
tex direction Qb for a given angle Q of the controlled 
variable H. To check whether this model describes the 
basic vortex lattice behavior when H is tilted away from 
the c-axis, we will use the CD as internal field detec- 
tors taking profit from the fact that Jc maximizes when 
Qb = Qd- Thus, if we know Qjj and Qmax (the angular 
position of the maximum in Jc), we are able to deter- 
mine the misalignment Qb — O = Od — Qmax- Although 
the sign of such misalignment is solely determined by 
the sign of /, its magnitude also depends on additional 
factors such as sin(20B) and k^. Besides that, the mis- 
alignment is strongly temperature and field dependent. 
It is easy to see from eq. (||) that 0b — > 8 for large 
enough h. Throughout this paper we will change each of 
these factors in order to show that eq. (^) satisfactorily 
accounts for the observed properties. 



At this point it is important to note that, although 
the misalignment between B and H is a low field effect, 
eq. can only be used in the field range where eq. is 
valid, i.e., for H ^ Hd- It turns out that all our data are 
well described by eq. (^. However, the very dilute vortex 
limit is conceptually interesting, and for completeness we 
will discuss it in the last section. 



III. EXPERIMENTAL DETAILS 

The Y:123 single crystals used in this workpsiere grown 
by the self-flux method as described in r ef.p and ex- 
hibit a critical temperature of Tc — 92K. We also per- 
formed measurements on single crystals of the layered su- 
perconductor NbSe2 with Tc = 7.2K. In all cases, colum- 
nar defects off the c-axis were introduced by irradiation 
with SOOMev Au'^^+ ions, using the TANDAR acceler- 
ator facility (Buenos Aires, Argentina). Table I sum- 
marizes the information about geometrical dimensions, 
mass anisotropy, dose-equivalent matching field i?$ and 
angle Qd the CD with respect to the c-axis, for all the 
crystals measured. 

The dc magnetization measurements were performed 
in a commercial SQUID magnetometer with two pickup 
coils, and both components (longitudinal Mi and trans- 
verse Mt) were recorded. Samples could be rotated 
around the axis perpendicular to both the c-axis and 
the CD. Two measuring procedures were used. The first 
one consisted of collecting a set of isothermal magneti- 
zation loops, each one recorded at a fixed angle £L In 
the second one, that we developed more recently,E3 the 
sample is rotated in steps at fixed T and H, and M is 
recorded at a dense set of orientations (the sample is 
not rotating during the measuring scan) . By appropriate 
proccessing, both methods allows us to obtain the irre- 
versible masHpt|i|aa,tion vector Mi(i/, T, 9), as described 
previouslya'll2Hl3Jl3 . 

In general, the relation (given by the critical state Bean 
model) between Mi and the various components of the 
anisotropic Jc is complicated. However, for thin plate- 



like samples, and as long as tan 9 < Ly/t 



(which 



in all the crystals used in our study is true for almost 
all field orientations, except in a very narrow angular 
range near 9 = 90°), due to purely geonietrical con- 
strains the following conditions are satisfiediEJ First, the 
pinning-related persistent currents flow essentially par- 
allel to the ab-planes, thus only one component of Jc is 
involved. Second, Mj is almost normal to the sample 
surface. Third, the geometrical factor that relates the 
modulus Mi with Jc is almost independent of 9. As a 
consequence, a measure of Mi{Q, H, T) is equivalent to a 
measure of the in-plane Jc(9, iJ, T). 



IV. EXPERIMENTAL RESULTS 
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A. anisotropy-dominated case 

Figure 1 shows the angular dependence of the irre- 
versible magnetization for Y:123 crystals A and B at 
T = 70K for several fields. The most evident feature 
in this figure is the asymmetry of Mi(Q) around = 
arising in the uniaxial pinning due to the CD. For sam- 
ple A at H > lOkOe, we observe the well known peak at 
the tracks angular position 9_d = 32° , and at lower fields 
the peak progressively shifts away from toward the 
c-axis. A more complete set of curves showing the shiit 
for this crystal at various T and H can be found in ref.H. 
A similar behavior is observed in crystal B, although the 
shift turns out to be smaller than in A. These two crys- 
tals have the same anisotropy and irradiation conditions, 
but different shapes. Thus, at the same T and H all fac- 
tors in eq. (^) are identical, except for /(i>, e). As seen 
in table I, the difference in demagnetizing factors results 
in a smaller /(i>, e) for sample B than for A. Hence, the 
misalignment in sample B is expected to be smaller, as 
indeed observed. 



B. compensated case 

A striking result predicted by eq. (||) is that the com- 
peting effects (anisotropy and geometry) could be exactly 
compensated if one were able to tune the demagnetizing 
factors and the anisotropy in order to get /(i>, e) = 0, a 
condition that is satisfied for 1 — i^z = (1 — i^y) e^. For 
the Y:123 single crystals used in this work, with e ^ 1/7, 
this requires extremely thin samples with a big area. Ta- 
ble I shows that crystals C and D almost exactly satisfy 
this compensating condition, as the absolute values of 
their /(/>, e) are, respectively, a factor of ~ 20 and ~ 10 
smaller than in A. Fig. 2 shows the angular dependence 
of Mi for these two crystals at T = 60K. For the sake 
of clarity only a reduced set of fields is shown, and some 
curves have been shifted vertically. Since the CD orien- 
tations 0u in samples C and D are different, in order 
to compare them we set the abscise as the relative angle 
O — 8d . In this figure we clearly observe that the peaks 
remain locked at the tracks direction even for the lowest 
fields, in complete agreement with the expectation. The 
same behavior was observed for other temperatures. 

C. geometry-dominated case 

So far, the four samples studied were Y:123 crystals 
with the same anisotropy but different geometries. On 
these samples we observed that the peak either shifts 
in the direction of the c-axis or almost does not devi- 
ate from the CD direction. As pointed out previously, 
this behavior arises from the strong anisotropy effect in 
this material. In order to change the sign of the devi- 
ation (i.e., a shift toward the ab- plane), we need to re- 



duce the anisotropy effects. (Table I shows that crystal 
D has / < 0, thus strictly speaking it is in the geometry- 
dominated case, but the shift is too small to be detected) . 

To that end we decided to measure NbSe2 single crys- 
tals, which have e ~ 1/3, making the anisotropy effect 
about 5 times smaller than in Y:123. Besides this, very 
large and thin NbSe2 crystals can be readily found, and 
they can be easily cut to obtain the desired shape. So 
we irradiated a rectangular crystal (labelled as sample 
E), such that it is in the extreme geometry-dominated 
case, with /(z^, e) <C 0. Fig. 3(a) shows Mi{Q) for this 
sample aX T — AAK for several H. At high fields we ob- 
serve a peak at the tracks' direction Qd — 27°. As field 
decreases the peak becomes broader and, in contrast to 
the Y:123 observed behavior, it progressively moves away 
from Qd toward the ab-plane, in agreement with a neg- 
ative /{v, e). 

The conclusive evidence that the function /(i>, e) dom- 
inates the behavior of the misalignment Qb — © comes 
from samples F and G, which are pieces of crystal E. 
These samples were obtained by cutting the sample E 
along a line parallel to its shortest side, in such a way that 
the demagnetizing factor remains unaltered, but Vy in- 
creases. In this way, the absolute value of the function 
f{v,e) was progressively reduced, i.e., we moved away 
from the "extreme geometry-dominated case" and ap- 
proached the "compensated case" (see Table I). Accord- 
ing to eq. (|^), the deviation of the maximum in Mi(8) 
for given H and T should become progressively smaller 
for crystals F and G as compared to crystal E. This 
is in fact observed, as demonstrated in Fig. 4, where 
©max — 6d for crystals E, F and G is plotted as a func- 
tion of h = H/Hc2{T,e). 

The misalignments for all the Y:123 crystals shown 
in Figs. (1) and (2) are also included in Fig. 4. Thus, 
this figure summarizes all the samples studied in the 
present work, at various temperatures and fields. The 
three possible low field behaviors are clearly visible: 
anisotropy-dominated (upward curvature), geometry- 
dominated (downward curvature) and compensated (al- 
most horizontal curves). It is worth to note that, in 
all the non-compensated cases and for both materials, 
the misalignment between B and H becomes relevant for 
fields below a certain characteristic field H 0.02Hc2- 

D. Quantitative test of the model 

The above results clearly demonstrate that the quali- 
tative differences in the low field behavior are controlled 
by the factor /(i>, e). We now want to verify whether the 
H and T dependence of the shift is well described by the 
model. According to eq. (||) these two variables appear 
only through the combination h = H/ Hc2iQB, T). Thus, 
\&D — Qmax\ should increasc not only with decreasing 
H at fixed T, as already seen in figs. 1 and 3(a), but 
also with decreasing T at fixed H, due to the increase 
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in Hc2{T). This second expectation is also verified, as 
shown in Fig. 3(b) where Mi{Q) for sample E was plotted 
at constant field H = O.WkOe for several temperatures. 

The equivalence between the variations in T and H 
is quantitatively verified in the main panel of figure 5, 
where sm{Q,nax ~ 0d)//(^', e) is plotted as a function of 
h for the two sets of data shown in Fig. 3(a) and (b). We 
observe quite a good scaling, thus confirming that h is 
the appropriate variable. The upper critical field values 
Hc2{kG) = 51.8 (1 — t) were taken from the liter ature,t£l 
thus the superposition of the two curves involves no free 
parameters. 

Finally, we analize the quantitative effect of the factor 
/(i>, e). This factor is a constant for a given sample, so it 
is the same for all the data in the main panel of fig. 5. In 
contrast, in the inset we show the same scaling procedure 
for the crystals E, F and G at T = iAK, so now /(i>, e) is 
different for each sample, while all the other parameters 
remain identical. We again obtain a good superposition 
of the data, although the scaling is poorer than in the 
main panel, probably due to the damage produced in the 
crystal after each cut process. 

The solid line in the main panel of Fig. 5 depicts 
the expectation of eq. (||), with 9_b = 6_d = 27° 
as experimentally determined from the location of the 
maximum at high fields, and a single fitting parameter 
K — 5.6. The same curve is shown in the inset. Although, 
the value n = 5.6 is smaller than the accepted valueE£l 
K 9, a similar discrepancy was reported by Zhukov et 
al. when studying the lock-in effectia Y:123 by both twin 
boundaries] jand columnar defects. lij Moreover, in previ- 
ous studiesu^Ej in Y:123 and Er:123 with CD, we had also 
found that the results are well described using values of 
the penetration depth A smaller than the accepted ones. 
We will discuss this issue in the next section. 



V. THE INFLUENCE OF THE CD 

It is important to keep in mind that the shift in Qmax 
at low h is not due to the CD. We are only using them as 
a passive tool to measure the vortex direction in the bulk 
of the samples, which is not easy to do by other methods. 
In fact, the pinning of the CD is not contained in cq.(|^), 
which arises from the minimization of a free energy, and 
thus describes a state of thermodynamic equilibrium. It 
is obvious, on the other hand, that the uniaxial pinning of 
these correlated structures is relevant and sJjipuld be in- 
cluded in the analysis. This is usually dondaa by adding 
to the free energy (||) a term Fpin that accounts for the 
correction to the vortex self-energy due to the CD, and 
then comparing the energy of alternative configurations. 

This additional contribution depends on the orienta- 
tion of the vortices, Fpin = Fpin (0b)i Sind it is always 
negative, reflecting the fact that, for Qb ©d, the CD 
promote the formation of staircase vortices whose self en- 
ergy is lower than that of a straight vortex at the same av- 



erage orientation. Fpin (©s) decreases as Qb approaches 
Qd due to the increase of the core trapped fraction, and 
it minimizes for that orientation, when the vortex cores 
are totally trapped into the tracks The key point in the 
context of the present study, however, is that the incor- 
poration of Fpin (0s) into the scenario does not modify 
the previous results, as we show below. 

Let's first consider that H is applied at the angle Q = 
Qmax such that, in the absence of pinning and according 
to eq. (^), the vortices would be at the angle 6b = Qd- 
If we now "turn on" Fpin (©s), the only effect will be to 
deepen the already existing minimum of the free energy 
at this orientation, without changing the angle. 

Let's now suppose that H is applied at an angle Q 
slightly smaller or slightly larger than Qmax- In the ab- 
sence of pinning vortices would respectively orient at an- 
gles Qb slightly smaller or slightly larger than Qd, ac- 
cording to eq. (H). The addition of the term Fpin (Qb) 
will now shift the vortices towards Qd, that is, a kind 
of effective angular attractive potencial towards the CD 
orientation will develop. In particular, for \Q — Qmax\ < 
(fL, the influence of Fpin {Qb) will be so strong that the 
system will minimize its free energy by orienting the vor- 
tices exactly along the CD. This is the well known lock- 
in effect, namely that the internal field remains locked 
to Qd over a finite range of 0,rSi|B<ji fL{h) is called the 
lock-in angle. In previous worksoll3 we have extensively 
studied this effect, that manifests in our measurements 
as a plateau in Mi{Q)) of width 2ipL- Note that the center 
of the plateau coincides with Qmax- Thus, although the 
relation Qb vs 8 will be modified by the CD, the angle 
Qynax, experimentally defined as the maximum in Mi{Q) 
or as the center of the plateau where necessary, will still 
be given by eq. (^. 

As a final comment it is worth to mention that, al- 
though oc l//i, similarly to the misalignment between 
B and H, the physical origin of both effects is totally 
different. Indeed, in NbSe2 we observe the shift but not 
the plateau associated to the lock-in, probably because 
the pinning of the CD is less effective in this material as 
compared to the Y:123 and Er:123 compounds. 

VI. THE VERY LOW FIELD LIMIT 

For completeness, we now consider the very low field 
limit, although it does not apply to our data. In this limit 
only nearest neighbours (NN) vortex-vortex interactions 
need to be taken into account, thus 

^.^^^[1 + ,(B)] (4) 

where Hd (Qb) ~ H'^cI^q^. 

In an isotropic material, the sum of the 6 equal con- 
tributions from the NN gives (7(B) — kr~-^/^e~^ , where 

k — J^^y^ 1; r — a/X and a — \J'^ is the vortex 
lattice parameter. In this limit r ;3> 1, thus g (B) ^ 1. In 
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the anisotropic case, g (B) must be modified to account 
for the angular dependence of A, and for the distortion of 
the triangular lattice. Regardless of the details, it is clear 
that in this case g (B) will still be exponentially small at 
low enough vortex densities, and the same is true for the 
derivatives (dg/dBi). Thus, replacement of (jj) into eq. 
(|l|) gives 



h: 



sin6_B 



cl" 



[1 + % (B)] 



H,, 



VyBy 



cos9s , H 

— -[i+vz m = - 



1 - u,, 



— VzBz 



(5) 



(6) 



where r]i (B), which account for NN interactions, are ex- 
ponentially small in the limit of -B — » 0. 

By taking the ratio of both components, we can calcu- 
late the relation between and 8s for the first vortex 
to penetrate. 



e tan 9b 



(g, - VyBy) (1 - (1 ~ 

{H, - V,B;) (1 - Vy) ^ (1 - Vy) 



tan 6 (7) 



So, in the very dilute 

limit the system is in the geometry-dominated case for 

(\ ^ Vy) > (1 — Vz) and in the anisotropy-dominated 
case for {1 — Vy) < {1 — v^), exactly the same result 
that we had found for intermediate fields. 

Eq. (HI) predicts that the misalignment grows indef- 
initely as /i — > 0, even allowing for &b and 9 to lay 
in different quadrants. As an extreme case, that equa- 
tion has no real solution when the absolute value of its 
RHS becomes larger than unity. These clearly unphysi- 
cal results are just a manifestation of the inapplicabil- 
ity of eq. (^) at very low fields. As an example, we 
can estimate the field at which the peak in Mi{Q) is 
predicted to coincide with the ab-plane (9 = 90°) for 
sample E. The result is w 0.0013, or H/Hd ~ 0.09, 
[that is H oi 250e in fig 3(a)], clearly below the lower 
limit of validity. (Note that, due to the strong demagne- 
tizing effects, vortex penetration becomes energetically 
favourable not at Hd but rather at a much smaller field 
H*^ - {l-Vz)Hci - 0.02ifci for sample E.) In con- 
trast, eq. indicates that for this sample at very low 
fields the peak should be observed at 9 w 88°, as marked 
with an arrow in fig. 3(a). 

Finally, we must note that the influence of the uniax- 
ial pinning of the CD turns progressively stronger as h 
decreases and the term Fpi„ becomes a significant frac- 
tion of the total free energy. At low enough fields the 
lock-in angle (pL oc 1/h covers most of the angular range 
(this is essentially what Klein et al.El called flux-flop in 
their early work) thus producing a very broad plateau 
that turns unapplicable our method to determine 9maa:- 



VII. CONCLUSION 

In summary, we performed a detailed study of the in- 
fiuence of anisotropy and sample geometry in the deter- 
mination of the vortex orientation in bulk superconduc- 
tors. We showed that these effects are relevant at low 
fields and become negligible at high fields. On top of 
that, we developed a model that correctly accounts for 
the sample shape, mass anisotropy, temperature and field 
dependencies of the misalignment between the applied 
field and the vortex direction. We demonstrated that the 
sign of the misalignment is solely determined by a func- 
tion that contains the combined effects of the anisotropy 
and the demagnetizing factors, and we presented ex- 
perimental examples of the three possible situations, 
namely anisotropy-dominated, geometry-dominated and 
compensated cases. We finally discussed the very low 
vortex density limit. 

An important practical consequence of these results is 
that studies of the pinning properties of tilted CD that 
are based solely on measurements at H ||CD, or on com- 
parison of this orientation with a few others, give valid in- 
formation at high fields, but are misleading at low fields: 
vortices are just not oriented in the right direction. To 
avoid this problem, a rather complete knowledge of the 
angular dependent response is required. 
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TABLE I. Irradiation and sliape specifications for all the crystals studied. 
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